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It is generally accepted that doped Mott insulators can be well characterized by the t-J model. In 
the t-J model, the electron fractionalization is dictated by the phase string effect. We found that in 
the underdoped regime, the antiferromagnetic and superconducting phases are dual: in the former, 
holons are confined while spinous are deconfined, and vice versa in the latter. These two phases 
are separated by a novel phase, the so-called Bose-insulating phase, where both holons and spinous 
are deconfined. A pair of Wilson loops was found to constitute a complete set of order parameters 
determining this zero-temperature phase diagram. The quantum phase transitions between these 
phases are suggested to be of non-Landau-Ginzburg- Wilson type. 
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Introduction. -The concept of fractionalization, find- 
ing its analog in quantum chromodynamics, is nowa- 
days a guiding principle of strongly correlated systems 
[1]. Specifically, a quasiparticle (like electron), at short 
spacetime scales, is effectively fractionalized into a few 
degrees of freedom (like spin and charge) which, at large 
scales, are "glued" by emergent gauge degrees of free- 
dom. As such, the gauge fields affect profoundly the 
formation of various novel quantum phases [2]. A fun- 
damental issue is the nature of the quantum phase tran- 
sition between these unconventional phases. In the con- 
ventional Landau-Ginzburg- Wilson (LGW) symmetry- 
breaking paradigm, different phases are characterized by 
the presence or absence of certain local order parameters. 
On the contrary, the fluctuation of emergent gauge field 
may contribute to low energy excitations of the system. 
In general, the transition between such phases cannot 
be described by the symmetry breaking paradigm, as is 
exemplified in Ref. [3]. 

Practically, an important playground of fractionaliza- 
tion is doped cuprates [TJ |2] . Despite the debate about 
the nature of the cuprate phase diagram, it is widely ac- 
cepted that the essential physics of high Tc superconduc- 
tivity is captured by a doped Mott insulator where the 
electron fractionalization into spinous and holons is seem- 
ingly inevitable [J [^ HI [S] . Indeed, there have been abun- 
dant evidence indicating that the electron fractionaliza- 
tion leads to new quantum phases p] and their transitions 
may not be understood in terms of symmetry breaking. 
This opinion has been reiterated very recently in Ref. [B] , 
where it was conjectured that the drastic change in the 
nature of quantum statistics-a direct manifestation of 
electron fractionalization-is at the root of the pseudogap 
phase found in the cuprates. 

Many fundamental issues arise thereby. How does the 
electron fractionalization turn an antiferromagnet into a 
superconductor upon doping? What does the phase dia- 
gram look like and what are the underlying order param- 



eters? Substantial efforts [7HTT] suggest that an exact 
non-perturbative result, the so-called phase string effect, 
discovered for the t-J model [7] is at the core of these 
issues. In this Letter, we present an analytical study of 
these open issues based on this exact result. We stress 
that our theory may be extended to other systems such as 
the Hubbard model on the honeycomb lattice currently 
undergoing intense investigations ^12] . 
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FIG. 1: (Color online) Zero temperature phase diagram of 
underdoped Mott insulators. The ball (in blue) and the arrow 
(in red) stand for the holon and spinon, respectively. The 
vortex, in red (blue) , surrounding a holon (spinon) arises from 
the spinon (holon) condensate. The wavy line stands for the 
confinement. 



Main results and qualitative discussions. -h\ essence, 
the phase string effect renders an electron "fractional- 
ized" into two topological objects, the spinon and the 
holon, each of which is bosonic and carries a tt-Aux '13]. 
A C/(l) gauge field A'^ (^p, radiated by the holons 
(spinous), interacts with spinons (holons) through min- 
imal coupling. (This is the so-called mutual statistical 
interaction which was also found in different contexts 
[12l[T4j). The macroscopic electric current (density), j'', 
is fully carried by holons and driven by both the external 
electric field E and the "electric field" E* resulting from 
A'l induced by spinons. Interestingly, E'' finds its origin 
analogous to that of Ohmic dissipation in type-II super- 
conductors: each spinon mimics a "magnetic vortex" sus- 
pending in holon fluids and, upon moving, generates an 
electric field antiparallel to j'', i.e., E" = — tt^CTsJ'', with 



as the spinon conductivity characterizing the mobihty of 
spinons. From the Ohm's law, i.e., crj^^j'' = — 7r^(Tsj'' + E, 
with afi the holon conductivity, we find that the (electric) 
resistivity, defined as cr'^^j^ = E, is 



<T-'=a-' 



TT^CT,, 



(1) 



A microscopic derivation of this composition rule will be 
given later. 

For the doping 6 sufhcicntly small, the antiferromag- 
netic (AF) phase is reached where the spinons are de- 
confined and form superfluids, i.e., CTs — >■ c». According 
to Eq. ([I]) this phase is insulating, giving a vanishing a. 
Indeed, single holon cannot appear in the excitation spec- 
trum. Rather, holons are excited in pair and are logarith- 
mically confined. These pairs are bound to the vortices of 
spinon superfluids and are thus immobile. Upon increas- 
ing S, the typical size of holon pairs becomes larger and 
larger and, as a result, the confinement becomes weaker 
and weaker. Eventually, a quantum critical point (QCP) 
is reached. Beyond this QCP holons become deconfined. 
For sufficiently large 5, a dual scenario applies. The 
holons are deconfined and form superfluids, i.e., a^ — >■ cx). 
In this phase, single spinon cannot be excited. Rather, 
spinons are excited in pair, logarithmically confined and 
bound to the vortices of holon superfluids, i.e., Ug = 0. 
(Consequently, the excitation spectrum is composed of 
integer spin excitations.) According to Eq. ([T]) this phase 
is superconducting (SC), i.e., a — >■ oo. Upon decreasing 
S, the confinement becomes weaker and weaker. Even- 
tually, another QCP is reached: for smaller S spinons 
are deconfined. Thus, an intermediate phase can appear 
where both spinon and holon vortices are condensed. In 
this phase, as,h -^ oo implies ct = 0. This brings us to 
the term, Bose insulating (BI) phase. 

The zero temperature phase diagram described above 
is summarized in Fig. [l] All three phases are char- 
acterized by confinement or deconfinement of spinons 
and holons, formally implemented by unconventional or- 
der parameters-a pair of Wilson loops W|''[C] with C 
a spacetime rectangle. Specifically, they are defined as 
the expectation values of e* *c m ^ that probes the in- 
teraction between a pair of test holons (spinons). W|[C] 
(W^'[C]) displays nonanalyticity at the holon (spinon) de- 
confinement QCP. The existence of these nonlocal "or- 
der parameters" suggests the non-LGW nature of the 
quantum phase transitions: in contrast to an LOW type 
scenario of order parameter competing, the AF and SC 
phases are intrinsically incompatible such that they can- 
not be turned into unless confinement or deconfinement 
occurs. In particular, they are generally separated by 
the BI phase, with the AF (SC)-BI transition as a holon 
(spinon) deconfinement QCP. We now turn to present 
some technical details. 

Lattice field theory. -FoiTaally, we start from the Hamil- 
tonian of the t-J model, which consists of the superex- 
change and hopping terms, describing the spin-flip and 



charge hopping process, respectively. An exact transfor- 
mation j2] keeping track of the phase string effect trans- 
forms the superexchange term into {a — 1,2) 



Hj^ 



-ia-A 



^^'^C.&l..e-'^=(^)5.^..,6. 



and the hopping term into 

with the coefficients omitted and h.c. being the Hermi- 
tian conjugate. Here, a = x,y denotes the unit vector in 
the a-direction. The spinon {h\^^hi„,a ^fji) ^nd holon 
{hi, hi) operators are both bosonic. The mutual statis- 
tics obeyed by spinons and holons is accounted for the 
following topological constraints satisfied by the gauge 
fields. For each loop C in the lattice plane. 



E.ec ^iJc ^ ^ E.i„.idec (^^t - blh) mod 27r, 
E.ec A^c ^ ^E.i-.- h^h mod 2^. 



(2) 



: inside C ' 



Here, Jc{x) = +1 (—1) for the link x -^ x+a {x+a — > x) 
on C, and is zero otherwise. 

In the underdoped regime, it suffices to invoke the 
mean field approximation [71 |B], with the Hamiltonian 
simplified to il = -E»(^.E.e— 4=«&I_.6^^ + 
i;,e'(^=«+^"«)/i^/i,+a) +h.c. {Js,th > 0). The exter- 
nal gauge field A^{i) generates an electromagnetic field 
and couples to the holon degree of freedom. The Hamil- 
tonian H and the constraints (pi) constitute our exact 
starting point. 

To proceed further, we consider the coherent state 
path integral representation of the partition function, 
Z — Tre"^^, where f3 {—^ oo) is the inverse tempera- 
ture. The technical challenge arises mainly from taking 
the topological constraints (l2| into account. Such a task 
was fulfilled in Ref. [9^ in the continuum limit. Here, we 
extend substantially the previous results to a more real- 
istic lattice model. We are able to show Z — J"^ e~ ^^ ^ 
with 



^ = Cs + Ch + V 



^l^^.A (^. _ 2^^^.) d, ^A'{ - 2^.jV^^) . (3) 
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Here, ^^f''^ are integer fields and £s,h is the spinon and 
holon Lagrangian, respectively. 



t .,w < + ^' -•^^^' 



Ch - h^ (D^o + A" - thD'^) h. 



Here and below, the Einstein's summation convention 
is implied and the summation over the indices /i, v, A 
includes both the imaginary time and spatial compo- 
nents. The on-site repulsive potential V softens the 



hard-core boson condition and depends on h^h and bj^ba-. 
We shall not further present its details to which the re- 
sults below are insensitive, e'^'^'*' is the totally antisym- 
metric tensor. Finally, the notation "Jf^" stands for 
the integral over the fields: b^ ,b,h^ ,h,A''''^ and the La- 
grange multipliers A"''', and the summation over ^^'^j 



-iA" 



e-°e-'"°, £)g = do - i^g, D\^„ = a{do - i^lg) with d^ 



the lattice derivative. With these preparations, the pair 
of order parameters are defined as 



Ws'''[C] = Z- 



-i:a^-^K'''j^) 



(4) 



where C is a spacetime rectangle with length T (R) in 
the imaginary time (spatial) direction and T ^ R. 

Differing from the prototypical field theory [9 , J^ 
keeps firm track of the compact nature of A^;'^, which 
affects profoundly the ground state properties, as will be 
shown below. In particular, the lattice mutual Chern- 
Simons term, namely the last term in Eq. ([3]) (Such 
a term was found previously in a study of Josephson 
junction arrays [H].), is periodic under a shift: A^'^ — > 
A^^f" + 2TTm'^^,,yV^'^ -^ ,jV^'''' + 27rm^''* with m"^^ G Z. 
Moreover, summing up .yVi^' ' enforces e^^'^d^^Jj'* to be 

TTITT with TO S Z. 

Superconducting p/iase. -Consider the case of dilute 
spin excitations where we may ignore the spinon field, 
i.e., 6^ = & = 0. Then, A''|/ip -I- V gives rise to the holon 
superfluid. More precisely, factorizing the holon field as 
h{x) = |/i|e*^(^^ and inserting it into Lh + V, we find that 
the fluctuation of \h\ is massive, whereas the Goldstone 
mode 0{x) is massless. Therefore, we ignore the terms 
associated with the spatial fluctuations of \h\ and obtain 



Ch = iW{doe-Af,- 

in the absence of A] 



zA'')-2i„|/i|2^„cos(d„0-A^) 
which is further simplified to 



i\hWdoe - Ag - lA") -f 4|/i|2[(d«0 - A%Y - 2]. (By the 
definition of H, a 27r-shift in daO is absorbed into A^^.) 

To calculate W^'^p], we separate A^ into the back- 
ground value and the fluctuation. The former leads to 
a uniform flux, 7r(5, at each plaquette and does not con- 
tribute to W^ [C] . Then, we introduce the unitary gauge 
so as to incorporate d^O into A" , and insert the simpli- 
fied expression of Ch into jSf . Integrating out the matter 
and A'' fields, we find 
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upon appropriate rescaling, where F/^^, — d^a^ — '^i/Q^u is 
the Maxwell tensor with a^ the fluctuating component of 
A^ , and ^/2th5 is the bare "charge". In the subsequent 
step, we integrate out a'* by using the Feynman gauge, 
which leads to important consequences. First of all, we 
find that at the ground state (where the external source 



■kJ^ is absent), two phase vortices of the holon super- 
fluid carrying opposite vorticity ■iz2'Ke^^^^d^,yVJ' are log- 
arithmically confined. Therefore, no free phase vortices 
exist with jV^ set to zero. Then, ttJ^ mimics an exter- 
nal dipole which may be produced by a pair of excited 
spinous with identical or opposite spin polarizations. In 
the latter case, a vortex with a vorticity of — 27r is excited 
from the background and bound to a spinon, forming a 
dipole. Taking these considerations into account, we find 



\nWm 



-irtfidTliiR, 



(6) 



which shows that the spinous are logarithmically con- 
fined. To calculate W|[C], we ignore a^. Integrating out 
the matter fields, we find the holon deconfinement. 
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-tHSJ:^AlA'^'+^J:.^AU^ 



1 



(T + R). 



(7) 



To further probe the SC long range order, we con- 
sider the response to a static small magnetic field. For 
this purpose we need to substitute A^ by A^ -I- A^ in 
Ch- If A^ ~ 2rmr (to, € Z), the partition function Z re- 
mains unchanged. For A'^ = (277i-|-l)7r, the additional tt- 
phase may be absorbed into A^ , leaving Ch (and thereby 
Z) unaffected. According to Eq. pi, we conclude that 
a single spin, with two possible polarization directions, 
is locally excited and nucleated at the magnetic vortex 
core [TU]. As such, the external magnetic field rmr flux 
is fully screened-a profound result of the integer field 
^''. If the flux value is not nnr, the magnetic field is 
excluded by the superconductor. Indeed, in this case, Z 
is merely determined by A'^ modulo tt (denoted as A^), 
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^2th<5;3 5;:icCos(d„e-A=) 



-t/^S/SEJAKq)!^ 



which justifies the Meissner effect. Here A'^{q) is the 
transverse part of the Fourier transform of A"^ . 

Antiferromagnetic phase.-We turn now to the case 
where holons are dilute and, likewise, set h^ = h = 0. 
Then, A'^dfe^^p -I- |6^p) -I- V^ gives rise to a two-component 



spinon superfluid, with |6^| 



y'n implying a 



magnetization in the transverse direction. Here n is 
the concentration of condensed spinous depending on 
d. Similar to the discussions on the SC phase, we 
factorize the two-component spinon field, (6l(x), 6^(a;)), 
as (|6^|, |6^|)e~*^'^^'' and insert it into Cs + V. Ignor- 
ing the spatial fluctuations of |6|, we obtain Cs — 
Ea K\^H<^doe - aA^ - zA^) - 2J, E. cos(d„0 - A^J]. 
^ ) Integrating out the matter and A'^ fields gives 



w![c] 



{.A"-} 






where Ff, 



d^Ajj — d^A^. 



Eq. (pi) has far-reaching consequences. First of all, by 
integrating out the A^ field, we find that at the ground 
state (tt J^ = 0) the phase vortices of the spinon super- 
fluid carrying opposite vorticity ■iz2ne^^'^d^jVj^ are loga- 
rithmically confined. Most importantly, ttJq mimics an 
external dipole produced by a pair of the holon and anti- 
holon. The latter is a — 7r-fluxoid, formed out of a — 27r 
phase vortex and a 7r-flux carried by the holon. Such a 
holon-anti-holon pair is logarithmically confined. 



In W-|[C] 



-27rJ.nTlni?. 



(9) 



In other words, a holon pair nucleates in a phase vor- 
tex of the spinon supcrfluid of vorticity — 27r, forming 
a "neutral" object. Furthermore, since the phase vor- 
tex is static, the pair of holons is spatially localized and 
the AF phase is insulating (see below for further expla- 
nations). It should be noticed that without the integer 
field ^1^ , such an insulating phase cannot be established. 
Instead, the SC phase is pushed all the way to 5 = Q [5]. 
Repeating the derivation of Eq. ([7]), we find the spinon 
deconfinement , 



\nW^[C] 



1 



8J^n 



{T + R). 



(10) 



Bose insulating phase.-The analytic results, namely 
Eqs. (|6]), ([7|, ([9| and ( 10 ), allow us to make an important 
observation. The asymptotic behavior, Eqs. ^ and (10 1, 
signal a critical concentration separating the spinon con- 
finement and deconfinement phases, at which WJ*[C] is 
nonanalytic in 5. Indeed, in the SC phase, the spinon ex- 
citations become progressively important as S decreases: 
they cause a renormalization of the bare "charge" in 
Eq. (pi) and eventually drive the system to the spinon 
deconfinement QCP where the "charge" vanishes. For 
smaller i5 the SC long range order disappears. Likewise, 
Eqs. (It]) and ^ signal another critical concentration, at 
which Wg [C] is nonanalytic, separating the holon confine- 
ment and deconfinement phases. Renormalizing the bare 
"charge" in Eq. ^ drives the system towards this QCP. 
For S larger than this critical value, the AF long range or- 
der disappears. Because these two renormalization mech- 
anisms are independent, these two QCPs are generally 
not identical, giving an intermediate phase where both 
the AF and SC long range orders vanish and both the 
holon and spinon are deconfined (condensed). 

This intermediate phase does not support dc electric 
transport. To prove this, we notice that the composition 
rule M is valid for the entire underdoped regime. Indeed, 
minimizing ^ gives 



SC, 



i.e 
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27r,^,-.''), 
a.E!^ and 



By definition, j^ 
j'^ = (Th{E^ + Ea). Noticing that the electric and holon 
currents are identical, we obtain Eq. (fTl) from these three 
relations. The striking structure of the composition rule 



is a consequence of i) that the two pieces involved are 
vortices, and ii) that they obey mutual statistical inter- 
action. (This fact has been established in a completely 
different context [S].) For the BI phase, CTs — > oo be- 
cause of spinon condensation, implying cr = 0. 

Crossover from Mott's law to activation /aiw. -Finally, 
we discuss a possible experimental observation of resis- 
tivity in the AF phase at sufficiently low but nonvanish- 
ing temperatures. In this regime, a holon-anti-holon pair 
bound to the vortex of spinon superfluids displays two- 
dimensional variable-range hopping conduction known as 
the Mott's law, In cr^^ ~ T'^/^ [16]. On the other hand, 
spinon superfluid is actually the so-called strongly en- 
tangled vortex phase [T7]. Carrying the arguments of 
Ref. [T7] to the present context, we find that the spin 
transport, in response to E'', displays an activation law, 
i.e., Incr^ ~ T~^. According to Eq. (|l|, the resistivity dis- 
plays a crossover from Mott's law to activation law upon 
decreasing temperatures, which may serve as a probe of 
holon (spinon) confinement (deconfinement) in the AF 
phase. 
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